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Spin-manipulated doping with magnetic (Ni) and non-magnetic (Mg) dopants constitutes the 
experimental attempts to obtain a singlet ground state system from the linear chain Heisenberg 
antiferromagnetic Cu-based d 9 spin-1/2 trimer compound Ca3Cua(P04)4 with doublet ground state. 
The present study is a first-principles based investigation of the effects of such doping on the spin- 
exchange mechanism and electronic structure of the parent compound. Site-selective doping with 
zero-spin dopants like Mg is proved to be more efficient than an integral spin dopant Ni in obtaining 
a spin-gap system with singlet ground state, as also observed in the experimental studies. Doping 
induced dimerized state is found to be the lowest in ground-state energy. Calculated spin exchange 
values along various possible paths resemble nicely with earlier experimental results. 

PACS numbers: 75.50.Ee, 71. 15. Mb, 75.10.Jm 



I. INTRODUCTION 

Basic concepts of spin-exchange mechanism and ex- 
cited state properties of low dimensional magnetic sys- 
tems can be useful to provide basic models for some com- 
plex phenomena like high-temperature superconductiv- 
ity in cuprates and arsenides. [TJ Low-dimensionality of 
magnetic exchange properties when associated with an 
energy-gap in the spin-excitation spectrum is even more 
interesting due to it's intrinsic two-level behaviour, [21 13] 
as demanded for utilization of such systems as a future 
quantum computing device material. [1] The simplest of 
such system, single-dimensional spin-1/2 Heisenberg an- 
tiferromagnetic (HAFM)dimer chain is well-known to be 
deprived of long-range magnetic order. [5J [B] It's zero-spin 
singlet ground state is separated from the integral-spin 
triplet excited state with an energy gap, known as spin- 
gap. Above a certain critical field(ft. c ), the Zeeman split- 
ted triplet excited states starts populating leading to a 
field-induced Bose-Einstein condensation and rendering 
the ground-state above h c to be the Bose condensate of 
low-energy magnons.[7] 

Inelastic Neutron scattering experiments and succes- 
sive analysis of exchange integrals have enabled to ob- 
tain an understanding of many such systems where the 
magnetic interactions within a geometric cluster of atoms 
are orders of magnitude larger than the inter-cluster 
interactions. [SI |9] The ground-state properties of such 
systems are very much sensitive of intra and inter clus- 
ter interactions. One such quasi- two-dimensional weakly 
coupled system is the triangle spin-1/2 trimer system 
La4Cu3Mo30i2, which can be well modelled by isolated 
Heisenberg spin-triangle model. [TU] Another geometri- 
cally frustrated isolated half-integral spin-trimer system 
is Vanadium Hallide, where the consecutive frustrated 
spins in an isolated triangular plaquette are aligned at 
120 degrees with respect to each other. [IT] The zero- 



temperature ground state magnetic structures of all such 
kind of systems are very much sensitive of the intra and 
inter-trimer interactions. 

In the present study, the parent system, 
Ca 3 Cu3(P0 4 )4 is a model spin-1/2 HAFM spin- 
trimer system with two symmetric site of Cu, where 
Cu2-Cul-Cu2 trimers form a linear chain along fe-axis 
with intra and inter-trimer exchange to be 126 K and 
3 K respectively. [5] [T3] The ground-state for this 
system will be a doublet with magnetic configuration 
Cu2 2 +(t)-Cul 2 +(|)-Cu2 2 +(t). In principle, if Cul 2+ 
can be replaced by a spin-one ion like Ni 2+ , the ground 
state of the system may be converted into a singlet one. 
The experimental observation by Pomjakushin et. al. 
|14)with powder inelastic neutron scattering analysis 
reveals that such selective doping does not occur in 
reality. Stable structure with Ni substituting Cul belong 
to a different space group C2/c with a doubled unit 
cell along c-axis, which differs from the original space 
group P2i/a for the parent compound Ca3Cu3(P0 4 ) 4 . 
Instead of Cul, Ni prefers to replace the end position 
Cu2 of the trimer chain. Hence, Ni does not satisfy 
the criteria of becoming a suitable dopant to induce a 
singlet ground-state in the substituted compound. In a 
successive study, Ghosh et. al. could successfully obtain 
a singlet ground state for this compound by selectively 
substituting Cu2 with a non-magnetic ion Mg only at 
one end position of the trimer chain. |15j 

These two studies have provided the motivation for 
the present investigation on Ca3Cu3(P04)4 system un- 
der spin-manipulated doping with Ni and Mg, where 
with a detailed electronic structure calculation, we aim 
to understand the effects of these dopants on the spin- 
exchange behaviour and ground-state. Experimentally, 
when a system is substitutionally doped with a particular 
dopant, there can be various possibilities of magnetic con- 
figuration. For this particular system, only a few of such 



FIG. 1: (colour online) (a): In a double-trimer array, the square planar and distorted square pyramidal coordination 
with O-ligands for Cul and Cu2 ions. Four corners of the Cul-03-04 Square are shared with PIO4 tetrahedra. (b): 
The double trimer Cu2-Cul-Cu2 array (discussed in text). Inter-trimer exchanges are mediated via 04 of PIO4 
tetrahetra. (c) Schematic view of the significant exchange couplings in the 6-axis trimer array. Jintra = J\ an d 
Jinter = J2 (see text for details) are denoted by solid and dotted lines. 



configurations will result into a singlet ground state spin- 
gap system. With the help of first-principles study, we 
have explained the modification of spin-exchange mecha- 
nism and exchanges along various possible pathways due 
to doping, which resembles nicely with the earlier exper- 
imental studies. fPfl ITS'] 

The paper is organized as follows. In section II, we ex- 
plain the details of the first principles calculation, where 
in various subsections, the structural details and elec- 
tronic structure investigations of Ni and Mg-doped sys- 
tems are presented. Exchange integrals along different 
superexchange paths are calculated and compared with 
experimental values. Partial density of state (PDOS) 
plots and the relevant band-structural details are also 
studied to understand the magnetic response of the 
doped systems. The last section summarizes the results 
with a conclusion. 



II. FIRST PRINCIPLES CALCULATION 

Ca3Cu3(P04)4, a suitable candidate to study spin- 
trimer antiferromagnetism crystallizes in the monoclinic 
structure with space group P2\/a. The system consists 



of two symmetric site for Cu, viz. Cul and Cu2 having 
square planar and distorted square pyramidal coordina- 
tion with O-ligands, as depicted in figure 1(a). The intra- 
trimer superexchange path between Cul and Cu2 (3.45 
A) is mediated via 04 consisting of the phosphate tetra- 
hedra centering PI. Inter-trimer superexchange paths 
Cul-03-Cu2, Cul-03-Pl-04-Cu2 and Cul-04-Pl-03- 
Cul are longer than Cul-Cu2 inter-trimer direct dis- 
tance(3.57 A), as can be seen from figure 1(b). 

We have investigated the electronic structure of 
this system with the help of linear muffin-tin orbitals 
(LMTO)[lSJ [T7] method in the atomic sphere approxi- 
mation (ASA), with combined correction terms, within 
the framework of LSDA + Hubbard U approach. The 
Hubbard U is added explicitly to the LSDA Hamilto- 
nian using the methodology prescribed by Czyzyk and 
Sawatzkv.fTB"] For the exchange correlation potential, 
we have used Perdew-Wang parametrization. |19j Self- 
consistency is achieved by performing calculations on a 
8 x 16 x 12 A;-mesh for the Brillouin zone integration. 
The on-site d-d Coulomb interaction and exchange inter- 
action parameters ( U, J) for the transition metal (TM) 
ions Ni and Cu are taken to be (4.5, 0.8) eV and (5, 1) 
eV respectively for all of our bulk calculations. However, 
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Magnetic Configuration 



Ground state details(LSDA+U) 




Nill(AFM) 
u Ni = 1.81 U B , u Cul = 0.8 u B , u Cu2 = 0.76 u B 

J mtra Ni-cu2 = 15 - 97 K, J 

S T =0, % singlet, doublet 



rttra cui-cuz = 62.10 K 




Nil2(AFM) 
u Ni = 1.81 u B , u Cu2 = 0.76 u. B 

j intra = 15 97 K 

J Ni-Cu2 R 

S T =0, singlet 




Ni21(AFM) 
U Ni = 1.75 U-b, |l Cul = 0.8 U B , U Cu2 = 0.76 U B 

jintra = 14 1 K J intra - 51 72K 

J Ni-Cul J - H - 1 *V J Cul-Cu2 3J - -/ * ,V 

S T =1, triplet 




Ni22(AFM) 
u. Ni = 1.76 u. B , u Cul = 0.8 u B 

J intra Ni-cui =14.05 K 
S T =3/2, quintet 



Table I: (colour online) Various substitutional possibilities of Ni-substitution and their corresponding ground state details 

resulting from a LSD A + U calculations. 



these parameters are optimized after comparing the the- 
oretically obtained exchange coupling constants with the 
experimental results. In the successive sections, this will 
be discussed in detail. 

To obtain a double-trimer array along &-axis, as seen 
in figure 1(b), we have constructed a 1 x 3 x 1 supercell 
to calculate the intra and inter-trimer exchanges. Realis- 
tically, in such an array, experimental substitution with 
Ni and Mg may lead to partial or complete substitution 
at various Cul and Cu2 sites. In the present study, we 
investigate all possible magnetic configurations resulting 
from substitution of Ni or Mg in the parent system. We 
denote these feasible configurations by Nijk, where j is 
the site position of Ni-doping and k is the number of such 
sites replaced with Ni per trimer. Thus, in case of Ni- 
substitution, there can be four substitutional magnetic 
configurations of the double-trimer array in Figure 1(b), 
viz., (a) Mil - substitution at a single Cul of one of the 
trimer, (b) Nil2 - substitution at Cul of both the trimers, 
(c) Ni21 - substitution at a single Cu2 of the trimer, (d) 
Ni22 - substitution at Cu2 of both the trimers. For Mg- 
substitution, the corresponding configurations are named 
as Mgll, Mgl2, Mg21 and Mg22 respectively. In the next 
subsections, we will describe the ground state possibili- 
ties, corresponding intra and inter-trimer exchanges and 



the related electronic structure implications of all such 
configurations. 

The model Hamiltonian for this particular system con- 
sists of two types of exchange couplings, viz. Jintra = Ji 
and J inter — Ji- Simplified Heisenberg Hamiltonian for 
this particular system can be written as: 

n 

H = J\ ^ S3j_i .(S3i_2 + S 3 j) + 
i=l 

71 

Ji /J(S3i_l.S3i-)-i + S3j.S3i+2) (1) 
i=l 

Here n is the number of trimers in an array. Schematic 
diagram of the trimer array is presented in figure 1(c), 
where different TM-spins are numbered according to the 
model Hamiltonian. 



A. Substitutional doping with Ni: resulting 

systems 

For the undoped parent compound, inter-trimer Cul- 
Cu2 exchange, as resulted from a LSDA + U calculation 
is ~ 92 K with the inter-trimer value to be ~ 5 K. For 
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FIG. 2: (colour online) The PDOS for Cul-3d, Cu2-3d, 

Ni-3d, 04-p, P-p states are plotted for (a) Nill, (b) 
Nil2, (c) Ni21 and (d) M22 resulting from a LSDA + U 
calculations. 



the Ni-substituted compound, among the four possibili- 
ties, experimental investigations by Pomjakushin et. al. 
[H] have discarded Nil2, since the composition with all 
Cul sites substituted with Ni crystallizes in a different 
space-group C2/c with a doubled unit cell along c-axis. 
The symmetry analysis presented at reference |14j implies 
that both Cu and Ni spins can be treated in a collinear 
arrangement. In addition, experimental observation in- 
dicates that Ni is seen to preferentially substitute Cu2 
at the end positions of the trimer. We have employed 
LSDA + U method for computing the total energy and 
exchange coupling constants, since the magnetic proper- 
ties of such system are totally driven by highly localized 
TM-3d states at Fermi-level. Computation of total en- 
ergy by LSDA + U calculation results into stabilization 
of AFM ground state with all adjacent spins in antipar- 
allel combination for all the four Ni-substituted configu- 
rations. Interestingly, Nil2 and Ni21 is found to be the 
highest and lowest among them in ground-state energy. 
It may be mentioned in passing that by LSDA calcula- 
tion, the lowest energy configuration is Nill with a fer- 
romagnetic (FM) spin- alignment. Table I pictorially rep- 
resents various above-mentioned configurations and their 
corresponding ground-state details resulting from LSDA 
+ U calculation. Exchange coupling constant along a 
specified path is computed by equating the ground-state 
energy difference of the anti-parallel and parallel spin 
combination along that path with the energy differences 
calculated from the model Hamiltonian (equation (1)) 
with relevant terms incoporating individual calculated 
spin-values for the magnetic ions as mentioned at the 
table. Total spin St of each such trimer configuration is 
specified in the same table. 

For Nill, St will be either or 1/2 corresponding 



to a singlet or doublet ground state with zero-spin Ni- 
substituted trimer and spin-1/2 original one. The inter- 
trimer exchanges are found to be very small in compari- 
son with the intra-trimer ones and are therefore prone to 
small fluctuations. Nil2, although experimentally unfea- 
sible to obtain, may have a possibility of St — singlet 
ground state. In fact, an integral spin Ni dopant was cho- 
sen by the experimentalists keeping this configuration in 
mind. For Ni21 and Ni22, St = 1 and 3/2 corresponding 
to triplet and quintet ground states respectively. The 
intra-trimer Ni-Cu exchange is much smaller than the 
Cu-Cu (intra) exchange. 

Experimental values for the Ni-Cu and Cu-Cu intra- 
trimer exchange was found to be 54.98 K and 9.86 K. 
Corresponding theoretical values, as computed for Ni21 
(lowest energy configuration), are 51.72 K and 14.1 K 
respectively. In addition, comparing exchanges for Nill 
and Ni21 implies that, Ni, while replacing Cu2 low- 
ers the total energy of the system more than a Cul- 
replacement. Corresponding exchange couplings for Ni- 
Cul and Ni-Cu2 are 14.1 K and 15.97 K respectively. As 
expected, inter-trimer Ni-Cu and Cu-Cu exchanges are 
much smaller in magnitude having values of ~ IK and 
4K respectively. To investigate the nature of magnetic 
exchange in more detail, we have plotted the partial den- 
sity of states (PDOS) for some selective ions for all four 
configurations in Figure 2(a) - (d). Since, effect of 04 is 
the most prominent among the O-ligands, 0-p PDOS is 
plotted only for 04. 

Available states to provide superexchange related hop- 
ping paths are determined by the square planar and 
square pyramidal coordination with O-ligands around 
Cul and Cu2. For all these four configurations, Cul- 
and Cu2-3d states at the up and down-spin channels have 
x 2 — y 2 orbital-character, the orbitals being confined in 
the plane of O-square. For Nill, up-spin empty states 
are mostly of Ni-3d x 2_ y 2 and Sd^-orbital characters, 
whereas for Ni21, those are of Sd 3z 2_ 1 and 3d x 2_ y 2 or- 
bital character. In general, Ni-3ci states are situated at a 
higher energy in comparison to the Cu-3rf states, leading 
to lowering of the value of intra-trimer Ni-Cu exchanges 
than Cu-Cu ones, as resulted from the theoretical calcu- 
lation. Although Ni-Cu and Cu-Cu distance are of the 
same order, lesser probability of lowering the total spin 
of the system for Ni-Cu hopping leads to a decrease of 
the corresponding exchange couplings. Comparing figure 
2(a) with 2(c), the small disparity of Ni-Cul and Ni-Cu2 
exchanges will also be evident due to the location of the 
Cul-3d x 2_ J( 2 states at slightly higher energy. Magnetic 
moments of Ni, Cul and Cu2, as calculated for Ni21 are 
1.75 hb, 0.8 fj,B and 0.76 \ib respectively, whereas the 
corresponding experimental values are 1.8 /j,b, 0.57 [Ib 
and 0.56 fis- As mentioned in reference [14] . reduction of 
Cu-spins in the experimental case may be due to the frus- 
tration in the Cu-moment and due to a combination of 
multiple possible configurations. Magnetic moments for 
other configurations are mentioned in table I. Hence, as 
also obtained from experimental studies, Ni-substitution 
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Magnetic Configuration 



Ground state details(LSDA+U) 






Mgll(AFM) 

ntra - CJ fJC If Jintra = -2 12 K 

Cul-Cu2 >V J Oi2-Mg-Oi2 * 

S T =1,1/2 triplet, doublet 

Mgl2(FM) 

M-Mg = M-B' U Cu2 = 0.77 u. B 

lintra . 1 ill/ 

J Cu2-Mg-Cu2 ~ "*■•« * 

S T =1, triplet 
Mg21(AFM) 

M-Mg = M B / Hcui = 0-81u B , U Cu2 = 0.77 U B 
c ul -c U 2 = 54.06 K,J^ Cul . Cul =8.87K 
S T = 0, singlet 
Mg22(AFM) 
M- Mg = u B , n Cul = 0.82u B 



J 



Cul-Cu 



! = 8.87 K 



S T = 0, singlet 



Table II: (colour online) Various substitutional possibilities of Mg-substitution and their corresponding ground state details 

resulting from a LSDA + U calculations. 
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FIG. 3: (colour online) The PDOS for Gul-3d, Cu2-3d, 
04-p, P-p states are plotted for (a) Mgll, (b) Mgl2, (c) 
Mg21 and (d) Mg22 resulting from a LSDA + U 
calculations. 
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FIG. 4: (colour online) The PDOS for Mg-3s and 2p 
states are plotted for (a) Mgll, (b) Mgl2, (c) Mg21 and 
(d) Mg22 resulting from a LSDA + U calculations. 



in place of Cul or Cu2 does not lead to a singlet AFM 
ground state. However, the exchange coupling constants 



are very much dependent on the values of U and J pa- 
rameters and therefore an optimization of the U and J 
parameters are necessary. Optimized value of U and J 
parameters are chosen from a comparison with the ex- 
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FIG. 5: (colour online) Band structure and the 
corresponding total DOS for Mg21 are plotted from a 
LSDA + U calculations. 
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FIG. 6: (colour online) Band structure and the 
corresponding total DOS for Mg22 are plotted from a 
LSDA + U calculations. 



perimental value of exchanges. 



B. Substitutional doping with Mg: resulting 
systems 

Motivated by the study of reference jUj, Ghosh et. al. 
have succeeded in obtaining a singlet ground state by 
doping one non-magnetic dopant like Mg in place of Cu. 
Similar to Ni-substitution, experimental observation con- 
firms the tendency of Mg to substitute the end-position 
(Cu2)of the trimer. In a similar way as in the previous 
subsection, we theoretically investigate all four different 
Mg-substitutional configurations Mgll, Mgl2, Mg21 and 
Mg22 as presented in the pictorial representations of ta- 
ble II with an LSDA + U calculation. LSDA calculation 
stabilizes Mgll as the lowest energy configuration with 
ferromagnetic (FM) spin-alignment. 

All these systems stabilize in an AFM ground-state 
except Mgl2, where a FM ground state is dominating. 
Among these configurations, only Mg21 and Mg22 may 
lead to a singlet ground state with the former one be- 
ing the lowest in ground-state energy. For Mg21, Mg 
substitution transforms the trimer chain into a dimer 
chain. Resulting Cul-Cu2 dimer fe-axis array leads to 
a reduction of the total spin of this spin-1/2 Cu2-Cul- 
Cu2 trimer system and turns it into a St = sys- 
tem. Strong intra-trimer (Cul-Cu2) and weak inter- 
trimer (Cul-Cul)interactions determine the ground state 
properties of Mg21. In case of Mg22, the Cul-spin-1/2 
&-axis array have AFM-coupled consecutive Cul atoms, 
resulting into a spin-singlet ground state. In this case, 
inter-trimer interaction is the only possible one. Thus 
with Mg-substitution, the spin-manipulation effects can 



be achieved in a better manner in comparison to Ni. Cal- 
culated Cul-Cu2 intra-trimer and Cul-Cul inter-trimer 
exchange coupling constants are both of AFM nature and 
estimated to be 54.06 K and 8.87 K respectively. Experi- 
mental value for the intra-trimer exchange are 55 K [15] . 

While investigating the LSDA + U PDOS figure, pre- 
sented in Figure 3(a)-(d), both for Mgll and Mgl2 in 
AFM and FM states, empty spin-up channel is of Cu2- 
3d x 2_ y 2 character. For Mgll, empty states at spin- 
down channel is mostly of CvLl-3d x 2_ y 2 character. For 
the lowest energy configuration Mg21 with AFM spin- 
arrangement, hopping to the up-spin Cu2 and down-spin 
Cul-3d x 2_j,2 states will be more prominent. For Mg22, 
only Gul-3d x 2_ y 2 states are available for hopping. While 
investigating the effects of Mg-3s and 2p-states, we have 
plotted the corresponding PDOS in Figure 4(a)-(d) to un- 
derstand the effects of Mg-s and p states on the superex- 
change. In case of Mgll and Mgl2, effects of Mg-s states 
on the filled localized states of Cu are more, which makes 
the Mg bands more localized. Although very small, the 
Mg moments for Mgll and Mgl2 are slightly larger than 
those for Mg21 and Mg22. Thus, Mg-s states have some 
role in mediating the superexchange for Mgll and Mgl2. 
In the last two cases, Mg-s bands are less localized. There 
is a substantial amount of hybridization of Mg-s states 
with the P-p states for Mgll, Mgl2 and Mg21. For Mg21 
and Mg22, most of the magnetic hopping is controlled by 
Cul-04-Cu2 path. Mg -s and p states have, therefore, a 
strong hybridization with O-p and P-p situated around 
Mg-04-Cul path. 

In Mg22, only possible superexchange interactions 
along inter-trimer Cul-Cul array localizes the 04-2p 
states. 01 and 06-2p states are highly hybridized with 
Mg-3s and 2p states. For Mg21, 01 and 06-2p states are 
less hybridized with Mg-3s and 2p states in comparison 
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to Mg22, as can be seen from figure 3(d) and 4(d). 

In figures 5 and 6 , we present the LSDA + U band- 
structure and the corresponding total DOS for up and 
down-spin channels of Mg21 and Mg22, for which sin- 
glet ground states are possible. Position of highly local- 
ized empty Cv2-3d x 2_ y 2 bands in the up-spin and Cul- 
3d x 2_ y 2 bands in the down-spin channels (as seen from 
Figure 5(a))for Mg21 are situated at lowest energy among 
the four configurations, causing the easiest lowering of 
energy by superexchange interactions among these four 
configurations. For Mg22 in an AFM configuration, both 
spin channels are having empty Cul-3ii a ;2_j ( 2 bands. The 
band-gaps for Mg21 and Mg22 are ~ 1.3 and 2 eV respec- 
tively. Magnetic moments of Cul and Cu2 ions for Mg21 
are 0.81 [1b, 0.77 (1b respectively. Theoretical analy- 
sis of energetics along with the DOS and band-structure 
enables to achieve a detailed understanding of the stabi- 
lization of a singlet ground state by Mg-substitution. 

It may be also mentioned that structural distortion 
has very little to do with the magnetic exchange prop- 
erties of such trimer system. In the initial case, Ni 2+ 
being almost of the same size as Cu 2+ , does not affect 
the Cu2-04-Cul angle much. For Ni21, Ni-04-Cul an- 
gle is almost the same as Cu2-04-Cul angle. For Mg21, 
Mg 2+ having a larger size than Cu 2+ , Mg-04-Cul angle 
is slightly larger than Cu2-04-Cul angle. However this 
may only affect the exchange between the end points of a 



trimer (as in case of Mgl2) and has very little effects on 
the dimer-chain exchange properties of the lowest energy 
configuration Mg21. 



III. CONCLUSION 

In this study, we have explored the doping induced 
spin-manipulation of the Cu-based d 9 spin-1/2 trimer 
system Ca 3 Cu 3 (P04)4 with the help of first principles 
LSDA + U calculation to obtain a detailed description 
of all possible substitutional configurations and the re- 
sulting ground states. The exchange coupling constants 
along various possible paths are also calculated and found 
to resemble nicely with experiment. The analysis of DOS 
and band-structure enables to understand the electronic 
mechanism behind obtaining a singlet ground state by 
substitution of Mg instead of Ni. 
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